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Formules de nouveau genre appliquées à la correction d'une orbite parabolique. 


I. Théorie. ° 
Par Th. Banachiewicz. 
Généralités sur les formules de nouveau genre. 


$ 1. Nous entendons par „formules de nouveau genre“ les formules où figurent, comme éléments, 
les „matrices“ ; les matrices qui, dans ces formules, sont préférables au point de vue des calculs numé- 
riques, et que nous appelons, à défaut d'un autre terme, ,cracoviens“, diffèrent d'ailleurs substantiel- 
lement des matrices de la théorie classique de Cayley. 

Le but principal des formules de nouveau genre est de réduire l'effort mental du calculateur 
et d'économiser son temps. On y parvient, en première ligne, par une schématisation unique et générale 
des transformations linéaires, l'emplacement d’une quantité dans le schéma indiquant déjà suffisamment, 
sans que le calculateur ait à partager son attention entre la formule dont il a à se servir et l'exécution 
du calcul, les opérations élémentaires auxquelles cette quantité doit être soumise. L'utilité des formules de 
nouveau genre provient encore de ce qu’elles permettent souvent d'éviter l'emploi de quantités auxiliaires 
à déterminer d’après les tables trigonométriques. En outre — et ce sont déjà des avantages d'une autre 
nature — la structure simple de ces formules facilite aussi l'exposition et la compréhension de divers 
chapitres de l'astronomie théorique. Les cracoviens particuliers „rotationnels“ p q r permettent même 
d'écrire immédiatement la solution dans une catégorie de problèmes de l'astronomie sphérique que l'on 
n'obtenait que par des raisonnements plus ou moins étendus et en recourant aux formules de la trigo- 
nométrie sphérique (voir Suppl. Internat. Nr. 3 de l'Annuaire de l'Obs. de Cracovie, pag. 60—61, article 
Libratione optico de Luna). 

Il paraît que, depuis l'invention des déterminants, c'est à dire depuis près d’un siècle, on n'a pas 
introduit dans les calculs astronomiques de nouveaux symboles généraux !); encore les déterminants eux- 
mêmes, très importants pour la théorie, ne s'emploient que rarement dans la pratique. Un certain scepti- 
cisme a priori de la part des astronomes envers ,les formules de nouveau genre“ serait ainsi parfaite- 
ment naturel; nous croyons cependant pouvoir attendre, sur base d'expérience, que ce sentiment se dissi- 
pera dès qu'on essayera de les appliquer aux calculs numériques. 

Pour obvier à des malentendus remarquons encore que les formules de nouveau genre ne sont 
nullement une transformation des formules anciennes, malgré que l'on retomberait souvent sur les for- 
mules connues en les réduisant à la forme ordinaire; de même, une solution des équations linéaires par 
les déterminants ne pouvait pas être regardée comme une transformation de la solution ordinaire. On se 
gardera en général de transformer les formules de nouveau genre en formules 
anciennes; on arriverait peut-être quelquefois de cette manière à des relations jusqu'à présent incon- 
nues, mais on perdrait en même temps les avantages précieux de la nouvelle solution. 

En parlant de calculs nous avons toujours en vue les calculs modernes, c'est à dire effectués, 
s'il y a lieu, avec des arithmomètres. 


§ 2. Voici un aperçu des définitions et des propriétés fondamentales des „cracoviens“ (par abbré- 
viation C-s). 
Pour que les C-s a et b 


di di as1| bis bn Ozi 
a = dig 9 A39 b es bis Da m (1) 
âi3 À) 2,3) bis Drg Dag 


soient égaux, il faut et il suffit que les termes ayant les mêmes indices soient égaux, c'est à dire 
An = bua (i = 1, 2, 3; k= 1, 2, 3). (2) 


On distingue, comme dans les déterminants, les lignes et les colonnes d'un C-n, mais les lignes 
et les colonnes peuvent être en nombres inégaux. 


') Le calcul vectoriel permet dď’écrire diverses relations d'une façon condensée, mais pour le calcul numé- 
rique effectif ses formules doivent être transformées en formules ordinaires. 


La somme c =a + b de deux C-s est un G-n ayant pour membre général la somme des 
C-s a et b, c'est à dire 
Cir = An | Due (3) 
La différentielle d'un C-n a est un C-n da, ayant pour membre général la différentielle 
du membre général, Cest à dire da. 
On multiplie les C-s a et b (ou, si l'on le préférerait, on les compose) en formant les termes 
du C-n — produit c — a . b d'après la formule 


Ca = ân bn + do by + ds Dig (4) 


C'est à dire: le C-n produit a pour terme général la somme des produits des termes correspon- 
dants (se trouvant dans la même ligne) de la i-ème colonne du premier G-n par les termes de la k-ième 
colonne du second C-n. On multiplie donc colonnes par colonnes"). La multiplication des 
C-s, qui est une opération fondamentale dans les applications, rappelle le procédé, si simple et si familier 
aux astronomes, de la formation, dans la théorie des moindres carrés, des équations normales d'après les 
équations de condition. 

On désigne par -/ (le C-n unitaire) 


[1 0 0 
r, — () 1 0) (5) 
lo o 1l 
On trouve immédiatement 
1 Ayo ds] 
JA = ân än dy (6) 
dzi dy ass | 


et l’on voit que les lignes de a sont les colonnes de :./. a, ct les colonnes de a sont les lignes de 7. a. 
On écrit simplement 
J-a = o/a (7) 
et lon prononce œa „Inverse a“. 
On trouve d'ailleurs 


a.o =a, (8) 


la multiplication d'un C-n par le C-n unitaire ne change pas donc le C-n multiplié. 
Si a.b =c, ona 


b . A = e7 © (9) 


et l'on voit que le changement de l'ordre des deux facteurs change le produit en son „inverse“. 
On définit le produit de trois G-s par la formule 


a.b.c—{(a.b).c, (10) 


de même a.b.c.d—{a.b.c).d cte. pour un plus grand nombre des facteurs. 
La formule (9) peut-être généralisée comme il suit. 
Si 
a.b.c.d.e= K, albors e.ofd.-7c.oJb .a = -7K. (11) 


Nous ne démontrerons pas ici cette formule, très importante dans les applications. 


On démontre facilement que la formule ('1) est vraie pour (n + 1) facteurs, si clle est vraie pour n facteurs. 
Or, pour n = 3, il suffit, pour vérifier la formule (11), d'écrire explicitement le terme général de a.b.c et de c. Jb .a. 
Cela étant, la formule (11) est vraie pour un nombre quelconque des facteurs. 


Il en suit 
a.b.c—a.(c.-,b), a.b.c.d=a.{d.,}e.-/b) cte. (12) 


(la règle d'association) 


a.(b.c)=a.c7c.b a.(b.c.d)—=a., /d. /c.b etc. (13) 


(la règle de dissociation) 


') Cayley (Coll. Mathem. Papers, Vol. IT, memoir 152) compose les matrices d'une autre façon (il multiplie 
„les lignes par les colonnes“); ce procédé est un peu plus commode pour la théorie, mais il paraît encombrant dans 
les calculs numériques. En composant les matrices d'après Cayley le calculateur aurait à rechercher tout le temps les 
quantités à multiplier, celles-ci n'appartenant plus à la même ligne et étant dispersées dans le tableau ; cela deman- 
derait un effort mental et une attention continuelle et serait sans doute une source féconde d'erreurs. C'est pourquoi 
la théorie de Cayley ne pouvait pas être ulilisée par les calculateurs et semble même n'avoir pas été remarquée par 
les astronomes. La règle de la composition des matrices, employée ici, change profondément les propriétés des pro- 
duits de ces symboles, c’est ce qui entraîne la nécessité d’un nom spécial pour eux. Remarquons encore que Cayley 
n'introduit pas la notion de la différentielle d’une matrice. 
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Toutes ces formules s'obtiennent de la même façon d'après (11). Par exemple, si K=a.b.c, 
alors, d'après (11), &/K = c. -jb . a = (c . -Jb) . a, d'après la définition (10). En changeant l'ordre des 
facteurs on trouve Kæa.(c.,-7/b) cest À dire la première des équations (12). Pour démontrer la pre- 
mière des équations (13) on pose K = a . (b .c), d’où, en vertu de (9), oJ K = (b . c) . a = b . c . a, d'après 
la définition du produit. Il en vient, en changeant l’ordre des facteurs K = a . o/c . b, C. Q. F. D. Il sera 
utile de retenir cette démonstration pour taciliter la transformation des produits des C-s. 


$ 3. On pose encore 


j1 0 | í cosn 0 sin n ) {cos n — sin n Q 
p (2) =,0 cosn —sinn: q (n) =} 0 1 O ; r(n) ="sinn pi 0 (14) 
i0 sinn  cosn \—sinn 0 cosn’ \ 0 1} 


et l’on a évidemment 


oJ (n) = p(- n), e/A(n) = q (— n), oJT (n) =r (—n). (15) 


Les trois C-s p, q et r peuvent être déduits d’un seul par un changement convenable des lignes 


et des colonnes. 

Leur importance dérive du théorème suivant. Si x y z sont les coordonnées d'un point fixe M et 
si elles deviennent x’ y” z après la rotation (en sens positif) du trièdre des axes de l'angle n autour de 
laxe x, ou autour de laxe y, ou autour de laxe z, on a les relations fondamentales 


i- lol. p (n) i= fol. q (n) h= lal. r (n). is 


Si l’on a trois points fixes M, (x gi 21), Ma (x Yə 2), M3 (3 Yz z3), leurs coordonnées devien- 
nent, après la rotation n autour de l'axe X 


x” x xs) w A Xa 
2, U Y Ys (=i h Uni p (n). (17) 
z z Zi Z2 Z3 


L'équation (17), avec deux équations analogues pour l'axe y et l'axe z (ou une seule), peut ser- 
vir à une déduction en bloc des formules de la trigonométrie sphérique. 


Formation des coefficients différentiels pour une orbite parabolique. 


§ 4. Dans le problème de la correction d’une orbite parabolique approchée, on cherche à déter- 

miner l’ensemble des coefficients des expressions pour cos Ô dx, dö en fonction des différentielles des élé- 
LT P | TR $ cos ô dx, dÈ 

ments; nous désignerons cet ensemble par A à DE a di 

Nous désignons par dp, dq, dr les rotations élémentaires autour des axes ,orbitaux“ !) x y z 
(voir, pour la définition de ces axes, Circ. de l'Obs de Crac. Nr. 16) qu'il convient d'appliquer au triè- 
dre des axes °). 

On cherche les valeurs numériques des coefficients du tableau 


dq lakt wt de dq | dr 


cos da COS Ò dò 
A EES aie IEN OANE E E u Dirata AEE 18 
sec = Ca dkT, a dp, dq, dr (18) 
dû 


| 
Pour des raisons de symétrie on peut y adjoindre encore lexpression de la différentielle de la 


distance géocentrique dp. 

Nous allons déterminer d’abord les accroissements dx dy dz des coordonnées héliocentriques en 
fonction de dx dy dp dq ar. 

D'après la Circ. de l'Obs. de Cracovie Nr 16, on a pour une comète 


| EEIE -PPN 
u= tuf o o, QJ (19) 


1) C'est M. Gadomski qui a suggéré ce terme. 

2) Les différentielles dp et dq ont la même signification que dans le Babnbestimmung de M. Bauschinger, mais 
nous désignons par dr, pour des raisons de symétrie, le dS de l'astronome allemand. Si l'on retient dq auprès de dq, sans 
craindre la confusion, on peut, pour la même raison, employer dans les mêmes formules dr et dr pour désigner deux 
différentielles différentes. D'ailleurs la différentielle dr du rayon vecteur ne figure pas dans nos formules. 


En différentiant cette relation on trouve 


ral dx\ (P. P, PÀ fx ,[P. P, PÀ 
i ao La, Q, QJ ii Egin Q, QJ (20) 


Pour déterminer les différentielles des éléments vectoriels P.....Q., remarquons que les neuf cosi- 
nus P,....R, peuvent être considérés comme coordonnées rectangulaires, dans le système équatorial, des 
trois points P, Q, R, situés sur les axes orbitaux x y z à une distance + 1 de l’origine, P, P, P, étant 
les coordonnées de P etc. D'autre part, la relation (19) n’est qu’un cas particulier, pour z2—0, de la 
formule générale 


shaja fée à ) n 
Ee i -| M R, R.J 


qui relie les coordonnées orbitales avec les coordonnées équatoriales. On obtiendra donc d'après la for- 
mule (21), appliquée successivement aux P, Q, R, les coordonnées cherchées P, + dP,, P, + daP,, P, + dP, etc., 
si l’on y substitue, au lieu de x y z, les coordonnées, rapportées à la position initiale des axes orbi- 
taux, de la position qu’occupe P etc. après les trois rotations des axes dp, dq, dr. 

Les coordonnées de P, Q, R dans la position définitive des axes orbitaux, rapportées à ces axes 
définitifs, formant un ensemble -/, ces coordonnées formeront dans le système orbital initial, d'après des 
raisonnements analogues à ceux qui conduisent à l'équation (17), un C-n 


o7. r (—dr).q(—dq) .p(-ap) =r (dr). q(— dq). p (— dp), (22) 


parce que la position initiale des axes orbitaux peut être obtenue de la position définitive par une inver- 
sion de l’ordre et de la direction des rotations. 
En mettant (22) au lieu du premier G-n dans le second membre de (21) il vient 


$ AR Fa fr: Ej Pa) 
Q, Re =r (dr). q (— dq). p-ap). i0: ©, Q?» (23) 

i g R LR. R R) 

(cor igés) (init: 110 


ou bien, si l’on peut considérer dp dq dr commme différentielles 


PSP, SPA (ECC 0) 


: Q: Q, ET Q: Q, H medr e ee 4 (24) 
R, R, R, R, R, RJ -dq ap +1 
corrigés) (initiaux) 


cest ce qu’on peut écrire aussi 
[Pa B, Bi) eA, | 0 — dr + dq) 


diQ. Q, Qi—=10 Q Qiiar 0 — ap (24 vis) 
LE KR] UE LADITE dp ol” 
d’où, en particulier 
a P, PI ý pd a Ei > q (25) 
y r 
40. 0, 0j REBL TEE 


REMARQUE. Le lecteur qui serait peu familiarisé avec les raisonnements concernant les C-s, pourra vérifier 
facilement la relation fondamentale (24 bis), équivalente à neuf équations, à l'aide de considérations cinématiques. La 
différentielle dP, par exemple, provient d’un déplacement + dr, parallèle à y, et d’un deplacement - dq, parallèle à z. Il en 
suit dP, = dr cos (yx) - dq cos (2x) = dr . Q, — dq . R,, c'est ce qui est conforme à la relation (24 bis). 


Le second terme du second membre de la relation (20) peut donc s'écrire, en ayant égard à (13): 


(xl | 0 +ar -an |P- P, P, -ydr | PME P. 
Qi = wed Q- Q, Q: 
(F1 |— ar K + ap! M R, R i x dq + ha k, kl 
et la formule (20) devient 
jax| |> — ydr | |P: P P.) 
dy; dy + xdr Q- 0, Q. (26) 
| —xdq + | Pa 7 


Il y a lieu de s'occuper maintenant des différentielles des coordonnées polaires gćocentriques 
en fonction de dx dy dz. Ces différentielles dọ, pọ cos x, pdè peuvent être considérées comme les accrois- 
sements des coordonnées rectangulaires dans un système d'axes qu’on obtient par la rotation æ du trièdre 


Xyz autour de laxe z, suivie d’une rotation —®ò autour de la position changée de l'axe y. Cela étant, on 
conclut immédiatement d'après les relations (16) deux fois appliquées 


| la + 
frere del = e j) .r(“).q(—0), 
ou bien, en tenant compte de (26) | z 
dọ lax — y ar | di 11 419 
p cos Ò dxi =: dy + xdr | ây) 4): al 3 (27) 
| pd? | — xdq + ydp Fe azı 43 
où lon a posé 
äio dun an| |2: Q- z| 
Pa d1 el) = q(—8).r(- x). Le a 4 (28) 
d30 431 432 z z z 
En omettant dọ, en divisant par £ et en posant 
[an æra) | 0 — sin A] cos æ sin x 0) (P. Q. R.) 
É 4 apas == pe Pa | M sing pa y as À A (29) 
Xj Lo Cos ò: p. , À 
on obtient !) 
[cos 5 dal dx — ydr &#, | 
| d ò | = dy + xdr | La el (30) 
EE pan La An 


La formule fondamentale (30), où les &; se déterminent d’après (29), s'applique à toutes les sections 
coniques. Ce ne sont que les expressions dx dy en fonction des accroissements des éléments à, e, T qui 
dépendent du genre de l'orbite à corriger. 


$ 5. Pour la parabole on obtient, en posant 


dx dy 
CT -= a f a m = tgo,  M=k (t=T): V24 °} (31) 
les expressions suivantes de différents b -s 
FPT ee ca AE 1 RARE LEPER [m° — 0,2 m°l (32) 
LT bo] EE- E El ly — 2q| Lol r |r: q + 1,2 m’| 


Si de est insensible, on posera b, = b}, = 0. 


Nous ne nous arrêterons pas à déduire les relations (32); elles peuvent être vérifiées facilement 
par la relation 


fax\ __fdr\ | cosv sinvì 
Ga + trdo! — sin v cos v’ 


en empruntant les expressions de dr et rdv en fonction de dq, dT, de aux traités d'astronomie théorique 
par ex. Klinkerfues-Buchholz, pg. 694 et 700, Watson, chapter Il etc. 
Ceci posé, on a définitivement 


Ce < E è i 52 oi (33) 
dq, dRT, de, dp, dq, dr 2 21, 23 À o| Es | 
On AO 10 CF y À ent ds 
En somme, pour calculer les coefficients différentiels cherchés, on a à effec- 
tuer les opérations suivantes: on détermine les by, d’après (32), les «x; d'après (29) 
et l'on fait la multiplication (33). 


$ 6. Nous avons supposé les neuf constantes vectorielles P.....hR, connues, parce qu'on les 
obtient au cours-même de la détermination arithmométrique d'une orbite (voir Circ. de l'Obs. de Crac. 
Nr. 16). Toutefois des cas peuvent se présenter, surtout tant que l'emploi des formules de nouveau genre 
ne sera pas général, dans lesquels on aura à calculer ces constantes d'apres 6}, à, &, €. La formule de 
nouveau genre pour ces quantités s'obtient immédiatement. En effet, les coordonnées de P, Q, R dans le 


1) Pour une autre forme de l'équation (30) v. notre article On a certain mathematical notion and its astrono- 
mical applications (Bull. de l’Acad. Polon. de Cracovie, 1923), où l'on a rapporté les différentielles dx dy dz aux axes 
X y z en supposant ces axes fixes, et non mobiles, comme ci-dessus. La relation remarquable (24) fut découverte 
après l'apparition de cet article, 
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système orbital étant 1, 0, 0; 0, 1, 0; 0, 0, 1, c’est à dire o7, et les axe; orbitaux pouvant être mis en 
coïncidence avec les axes équatoriaux xyz par la rotation —® autour de z, ensuite —i autour de x, 
— 6} autour de z, et finalement —Ee autour de x, il s’en suit 


[2 Q 4 
2, Q, R = 2. r (0). p (i) .r(- 9). p-e), 
IP O R] 
ou 
P, Q 7 a 
12, Q, Rir = r (©). p (À . r (19) . p (—8) (34) 
LP Ọ à fe | 
Pour la vérification partielle on a la formule 
Eo pai | 
$ Q, "R 


le symbole 2 servant à indiquer la multiplication du G-n par soi-même. 


REMARQUE. ll sera peut-être instruclif de comparer la relation (34) à la solution ordinaire arithmométrique 
du problème de la détermination de P,.... Rz. 

On trouve — non plus par un raisonnement pur, comme l’on a obtenu (34), maïs par des calculs assez longs 
bien que faciles (voir pour les six cosinus P2... Q: le Journ. of the Brit. Astronom. Association, Vol. XXXII, p. 232) : 


Pz = cos w cos Q — sin w cos isin Q 

P, = (sin w cos i cos R + cos w sin R)cose — sin w sin į sin € 

P, = sin wsin į ccs £ + (sin w cos į cos R + cos w sin Q) sin e 

Q: = — (sin w cos R + cos w cos į sin £2) 

Q, = (cos w cos i cos & — sin w sin &) cos € — cos w sin į sin € (36) 
Q: = cos w sin į cos £ + (cos w cos į cos R — sin w sin QR)sin 

R: = sinisin Q 

R, = — sin į cos Q cos £ — cosisin s 

R: = — sin i cos à sine + cos i cos e 


Dans la formule de nouveau genre (34), déjà après l'introduclion des expressions explicites de p et r, il y a 
55 symboles significatifs ; dans les formules anciennes (36) il y en a 154. Le rapport 154: 55 = 2.8. 
Les équations de vérification, sous forme ancienne, seront 


P? + P’ +P} =1 P,Q, + P,Q, + P,Q, =0 
Q+Q+Q*=1 P,R,+P,R,+P,R,=0 (37) 
a i a QR,+Q,R,+QR; —0, 


contenant 60 symboles, contre 13 symboles de la relation équivalente (35), ou bien 22 symboles de la même relation, 
si l’on y écrit explicitement ‘7 


§ 7. On sait’) que Safford, Tietjen et d'autres savants ont réussi à simplifier les équations 
de condition par l'introduction de coordonnées polaires géocentriques rapportées à l'orbite (la „zweite 
Methode“ du Traité de M. Bauschinger). Peut-on obtenir les mêmes avantages en se servant de for- 
mules de nouveau genre ? 

Pour répondre à cette question on n'a qu’à multiplier les deux membres de l'équation (30) par 

4328 As `s) 


Í ` L » 9 2 
| š où l'on a posé s= Va + a. 
ieai Gris” j Ee 


On posera 
fcos À dh\ _fcosôdal f a&:s a: 
4 (38) 
lodhi PI a Re) 
[Pr Pra| b. aa S PA :s| farn ai al (39) 
Nm" ES : : 
lo s | |- az -S yo : s| |æ X2 Azz 
et la relation (33) devient 
cos X db, dX pr 12 = 2 0 AC lis 7 
dq dkT de dp dq dr [On 2z ds £ dd ui 
$ i Lo 0 0 +y —x à" Lo s | 


On aura donc à calculer, comme auparavant, les b, d'après les relations (32), les &, par la rela- 
tion (29), après quoi les relations (38) et (39) donneront cos X dp, d7. et Ba et l'équation (40) fournira 
les dérivées cherchées. L’équation en cos X db ne contiendra que 4 inconnues dq, dkT, de, dr, et les 
coefficients de ces inconnues dans l'expression de d% seront de l’ordre de grandeur de sin %, c'est à dire 


1) Voir Radau, Réflexions sur les métbodes de correction des orbites, Bull. Astron., 1904. 


petits pour les orbites peu inclinées à l'écliptique. De cette façon on obtient les mêmes résultats que 
ceux que donne la méthode Safford-Tietjen, mais on évite le calcul des coordonnées X et Ÿ et de 
l'angle parallactique W, et les équations ont la même forme que dans la méthode usuelle. 


$ 9. Après la détermination de dp dq dr on aura les éléments vectoriels corrigés par la relation 
(24). Ces éléments vectoriels jouant aussi le rôle des constantes de Gauss, le calcul, d'après la formule 
(19), de l'éphéméride corrigée, est ainsi grandement simplifié en comparaison avec les méthodes jusqu'à 
présent en usage. Il y a lieu de remarquer que dans la formule (24) les éléments vectoriels peuvent être 
rapportés aussi bien à l'équateur qu'à l'écliptique. 

Si l'on désire calculer les éléments traditionnels ©, i, © corrigés, on peut le faire de deux ma- 
nières différentes. Ou bien on les obtiendra d'après les éléments vectoriels nouveaux, en les tirant des 


formules 
R,=— sin Ísin Q, R,= — sinicos Q, R 
aIr A or 4 kR,=—= cosi, (41) 
P, = sin i sin 6, Q, =  sinicos &, 
ou bien on emploiera les équations bien connues 
di (dp cos & sin® — sin ð) 
PE OS — (A | — sin © COS ®&  — cos Hi (42) 
tgi dð dr 0 0 tg i 


Les relations (24) et (42) supposent dp dq dr infinitésimaux. Toutefois la seconde méthode (42) 
est moins générale, parce qu'elle demande que di} et d& soient aussi infinitésimaux. 

Si l'on ne pouvait pas négliger les carrés et les produits des «p dq dr, leurs valeurs, détermi- 
nées d'après les équations de condition, ne seraient pas elles-mêmes définitives. On prendrait alors pour 
point de départ des calculs nouveaux les éléments vectoriels d'après la formule (23) qui donne les cosi- 
nus satisfaisant aux conditions géométriques (35). 


Les avantages de la nouvelle solution. 


Nous allons résumer ici brièvement les avantages de l'application des formules de nouveau genre 
au problème qui nous occupe. 


1) Les relations (32), (29), (33) constituant l'ensemble des formules dont on a besoin, sont bien 
plus simples et plus claires que les équations ordinaires. Le nombre des symboles mathématiques figu- 
rant dans la nouvelle solution est 140 contre 376 symboles de la solution du Traité de M. Bauschin- 
ger. Rapport 376 : 140 = 2.7, en faveur de la nouvelle solution. 


Ces chiffres ne sont d’ailleurs qu’approximatifs, vu une certaine indétermination de ce qu'on doit considérer 
comme un Symbole. Les détails de notre compte sont: relation (29) 42 symb., (32) 54 symb., (33) 44 symb., total 
140 symboles. Nous n'avons pas compté les zéros, mais d'aulre part nous comptons les «,, qu’on pourrait omeltre tous. 
Dans le Bahnbestimmung der Himmelskôrper de M. Bauschinger lesb Bc C (pg. 455) contiennent 86 sym- 
boles, j J h, H (pg. 460) 66 symboles, cos à dx 113 symb., dè 111 symboles (pg. 460), total 376 symboles. 


2) Le calcul d'après les formules nouvelles ne demande d'attention que dans la détermination 
des valeurs numériques des termes des cracoviens; ces termes sont d'ailleurs en général bien simples. 
Une fois ces termes écrits, le reste du calcul est complétement uniforme et mécanique, et ne demande 
que 20 minutes pour la détermination avec 4 chiffres des 12 coefficients cherchés pour une position 
complète, en supposant un calculateur non entraîné et un arithmomètre non électrique. On évite complé- 
tement les angles auxiliaires. En somme, l'effort mental est minime et le temps nécessaire bien réduit. 


3) On passe de la solution en coordonnées habituelles à la solution d'après SaffordTietjen 
par une simple multiplication par un cracovien. 


4) On n’a pas à déterminer les constantes de Gauss qui disparaissent entièrement des calculs ; 
une fois qu'on a corrigé l'orbite, on dispose des constantes nécessaires pour le calcul d'une nouvelle 
éphéméride. 

5) Les formules sont arithmométriques. 

Quant aux éléments vectoriels figurant dans la solution, rappelons qu'on les obtient déjà au cours 


d'une détermination rationnelle de l'orbite approchée et que 6 éléments interviennent dans le calcul 
arithmométrique d'une éphéméride. Ces éléments sont des simples fonctions des constantes de Gauss 


Cracovie, 1925 Mai 13. Th. Banachiewicz. 


lII. Calculs. 
Par J. Witkowski et K. Kordylewski 
Nous avons pris pour point de départ de nos calculs la moyenne arithmétique des deux systè- 


mes d'éléments de la Comète Orkisz (1925C) de Möller et Johannsen (Circ. No 66 de l'Obser- 
votoire de Copenhague) et de G. Merton (Circ. No 5 of British Astr. Assoc.) 


TASo eiaa P, sQ +0:532148 — 0532847 —0:657943 . 
i P, Q, R,į = ;—0-796185 —0050665 —0:602924 1925.0 
q = 1.109541 |P Q, R | -+0:287930 --0-844602 — 0451206] 


Le calcul est basé sur les observations de la Comète depuis 1925 Avril 5 jusqu'au Mai 27, pu- 
bliées dans les À. N. 224 12, 15 et 16, Circ. d. Bur. Centr. Astron. à Copenhague No. 63, 65, 66, Beob. 
Zirk. 16, 17, 18, 19, Harv. Bull. 817, Nature 2898, ainsi que sur les observations, pas encore publićes, 
faites aux observatoires de Cracovie (23 obs.), de Varsovie (6 obs.) et de Poznan (6 obs.). 

Après lomission de 6 observations évidemment erronées, on a pu utiliser pour le calcul 80 ob- 
servations, qui dûment réduites pour la parallaxe etc. ont été comparées avec une éphéméride calculée 
avec 5 decimales. 

Les équations de condition ont été calculćes d’après les formules ingénieuses, exposćes dans la 
I. partie de cette Circulaire. Voici un exemple de ce calcul pour 1925 Avril 7:0 

x = 22h 28m 37 3 = +18036.0 p = 1.6942 M = 0.05748 m = +0.0574 r = 1.1132 


2; | 0 — 1883] | +9215 — 3883 0| | +5321 —5328 — 6579 [3111 + 106| 
as, Qat = < +5903 OÙ J 43883 +9215 OÙ 17962 — 507 —6029! = 11497 +5613! (29) 
leu asa) alaian 0 5505 iaia No 0 10000] |+2879 +8447 —4512 — 6284 —1324 
fbi, bel _ f+9013 O1 NS + 12741  f+10066 + 768] fbis\ — 317 1 0002 OU ig 
és def AF 773 +6030) F 1274 22100 ET. 15381 TA \10033/ — i8) P 
+1C066 + 768 0 0 Ou 1274) [3111 + 106 dq d4T de «dp dq dr 
— 567 —13381 +318 0 O +1109% 4—1497 +5613; = f- 3026 +1764 — 48 — 610 +5294 —1259) (33 
| 0 0 O +1274 —11059 oj |—6284 —1824 \— 212 —7503 +178 —232 +2017 +6194; ©?) 
Les suivantes équations de condition ont été obtenues de cette manière. 
O— C v Date n 


équations en z% 


dq — 303 —313 - 331 —382 —398 —452 — 568 — 757 — 3.5) —1.7 Avr. 70 18 
dkT [+17 +166 +142 + 85 + 68 + 12 —105 -295 |- 1] — 1.2 9.0 13 
de _— 5— 6 10— 16—= 180 25010 + 0.51 +0.7 13.0 10 
up — 61 — 84 —131 —226 — 249 —314 — 390 —355 TU + 1.4 210 11 

| dq | +529 +534 +540 +540 +537 +518 +450 +294 +140!  +0.2 23.0 10 
dr —126 — 109 — 72 + 18 + 414 +131 +310 +585 91.91 +13 29.0 9 
+39.6! —0.7 Mai. 9.0 7 

+542] +04 20.0 4 

équations en à 

dq '— 21 — 52 — 116 —255 —291 —393 —511 —100 — 341 —0.5 Avr. 70 13 
dk T — 750 — 760 —798 —847 — 856 — 866 —821 = 51) |- 0.6} —0.6 90 13 
de + 18 + 25 + 40 + 70+ 78 + 98 +123 +129 + 401 +0.2 13.0 10 
(dp — 23 — 30 — 39 — 35— 28 + 10 +140 4-426 jengo — 1.2 210 11 

| dq | [20 +190 +162 + 85 + 61 — 17 —162 -338| +13.0} +4.2 23.0 9 
dr +619 +634 +665 + 729 +744 +782 +796 +613 +2271 +0.7 29.0 9 
+382.7 —6.1 Mai. 9.0 7 

+ 3.4) +35 21.5 5 


En résolvant ces équations d'après la méthode des moindres carrés et leurs assignant un poids 
égal, on trouve d'abord de = + 0.00048 + 0.00092. En posant donc de = 0 on a pour les autres quantités: 


dy = —"000243 : dT = +™000802 ; dp = —'000361 ; dq = — 000224 ; dr = +-000038 
Les corrections des elements vectoriels ont été calculées à l’aide de la formule differentielle : 
| +0.532148 —0.796185 +-0.287 :30] | 1000000 — 000038 — 000224 +0.531980 — 0.796322 +0.287861 
— 0.532847 0.030663 +0.844692! ++ 000038 1000000 + 000361! = 10532630 —0.050417 +0.844844! (21) 
| 0.657943 0.602924 —0.451206] |+ 000224 — 000361 1000000 — 0.658255 0.602764 —0450066| 
Nous trouvons pour les éléments corrigés : 
= P, Q, R 0.53198 —0.53263 —0.65826] 
— 1925 Avril 1.47820 z + M +0. 53263 —0.65 
$ 6 l Q REL = 110.79632? —0.05042 —06076! 1925.0 
q = 1.10930 PQR +0.28786 0.844814 — 0.45097 


à 
Un calcul de contrôle effectué pour Avril 7, 23 et Mai 9 prouve que les residus provenant des 

éléments corrigés ne different des résidus des équations de condition que dans les limites des erreurs du 

calcul à cinq décimales. La valeur moyenne absolue de la difference des residus en question est 1”. 
Les éléments eorrigés ont servi à calculer l'éphéméride suivante : 


1925 Juin 40 du = 615mb z p = +81038.F r = 14982 p = 1.7213 
Juillet 6.0 10 10 15 +60 57.3 1.8451 2.8139 
Août 7.0 10 54 24 +47 48.9 22123 2.9167 
Septembre 8.0 11 23 59 +40 10.6 2.5806 3.3471 


Contenu du Nr. 17: L'orbite corrigée de la Comète Orkisz (1925 C). Formules de nouveau genre appliquées à la correction 
d'une orbite parabolique. 1. Théorie (par Th. Banachiewicz). WU. Calculs (par J. Witkowski et K. Kordylewski). 


Drukowano na zlecenie Obserwatorjum Krakowskiego. 
Dyrektor: Tad. Banachiewicz. 


